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Neon  ion  embedded  in  a  hydrogen  plasma  at 
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Energy-level  spectrums  for  a  Neon  ion  in  a 
hydrogen  plasma. 


SECTION  1 


INTRODUCTION 

Equation  of  state  data  (pressure,  total  energy, 
specific  heat  and  other  thermodynamic  properties)  for 
plasmas  in  a  wide  range  of  conditions  is  a  problem  area 
of  considerable  importance.  For  a  low  density  fully 
ionized  plasma  at  very  high  temperatures,  classical  EOS- 
results  are  available  and  are  generally  valid.  However, 
EOS-properties  for  high  density,  low  temperature  plasmas 
(including  the  "cold"  compressed  matter  phase)  are  much 
more  complicated,  particularly  for  partially  ionized 
conditions.  In  this  domain,  available  EOS-data  is  scanty, 
and,  in  cases  where  it  is  available,  the  models^  used  do 
not  adequately  take  into  account  the  physics  of  the  inter¬ 
acting  plasmas.  This  is  particularly  true  for  the  elec¬ 
tronic  contribution  to  the  EOS.  Depending  on  the  density 
and  temperature,  quantum  effects  are  very  significant  for 
electron-electron  and  electron-ion  interactions.  At 
intermediate  densities  and  temperatures,  plasmas  can  be 
"strongly  coupled"  as  well  as  partially  ionized.  Both  of 
these  aspects  need  careful  consideration  and  have  not  been 
properly  included  in  previous  models.  In  addition,  for 
improved  accuracy,  the  "cold"  compressed  matter  phase  need 
consideration  of  solid-state  effects. 

Plasmas  in  the  above-mentioned  conditions  occur  in  a 
variety  of  experimental  conditions--shock-compressed 
materials,  the  superdense  core  of  laser-fusion  plasmas, 
laser-imploded  solid  targets  or  multi-layered  materials, 
etc.  The  ongoing  research  on  these  processes  requires  a 
knowledge  of  the  appropriate  EOS-properties.  Improved  EOS 
data  is  also  needed  for  more  realistic  hydrodynamical 
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calculations  (including  radiation  transport) ,  for  study¬ 
ing  implosion  characteristics  of  cylindrical  radiating 
plasmas,  as  well  as  for  a  variety  of  atomic  calculations. 

This  need  motivated  a  more  realistic  EOS-model 

suitable  for  a  wide  range  of  plasma  conditions,  densities 

and  temperatures.  The  present  report  provides  such  a 

comprehensive  EOS-model,  justifies  the  theoretical  basis 

and  presents  numerical  results  obtained.  This  model  is 

2 

based  on  the  density  functional  method  generalized  to 
arbitrary  temperatures.  This  method  enables  a  systematic 
incorporation  of  the  above-mentioned  effects.  In  Section 
2,  physical  considerations  for  this  model  are  outlined, 
its  theoretical  basis  is  justified  and  the  implications 
of  the  numerical  results  are  discussed.  The  specific 
mathematical  details  of  the  scheme  are  given  in  Appendix 
I. 

In  Section  3,  research  on  problems  involving  polari¬ 
zation  shifts  in  dense,  hot  plasmas  is  presented.  This 
is  an  important  effect  that  affects  atomic  processes  in 
plasmas,  particularly  the  radiation  characteristics. 

These  atomic  processes  involve  consideration  of  the  number 
and  spacing  of  ionic  energy  levels,  level  populations, 
transition  probabilities,  etc.  These  parameters  as  well 
as  radiation  transport  through  hot  plasmas  are  known  to 
be  drastically  affected  by  a  dense  plasma  environment. 

The  plasma  polarization  shift  (PPS)  plays  a  significant 
role  in  this  context.  The  dense  plasma  environment 
produces  significant  lowering  of  the  ionization  potential 
which  can  eliminate  excited  atomic  energy  levels  and 
possibly  entire  stages  of  ionization.  Although  previous 
work  has  been  done  to  include  this  "lowering"  in  atomic 
calculations''",  the  majority  of  these  models  are  semi- 
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empirical.  Detailed  investigation  of  the  PPS  over  a  wide 
range  of  plasma  densities  and  temperatures  is  important 
for  understanding  the  radiation  characteristics  of  hot 
plasmas  as  well  as  for  a  variety  of  atomic  calculations. 

A  systematic  study  of  plasma  polarization  effects 
is  outlined  in  Section  3.  Our  model  is  discussed  and 
compared  with  other  commonly  used  static  screening  models. 
The  improvements  of  the  present  model  are  that  it  treats 
nonlinear  screening  effects  exactly  and  quantum  many-body 
effects  are  included  via  a  self-consistent  scheme.  This 
model  is  applied  to  Neon  and  Argon  impurities  in  dense, 
hot  plasmas  and  results  are  compared  with  those  from  the 
Debye-Huckel  model.  Drastic  differences  in  the  energy- 
level  structure  indicate  the  usefulness  of  this  nonlinear 
self-consistent  model  for  dense  plasmas  where  Debye- 
Huckel-type  models  are  inadequate. 

For  purposes  of  plasma  diagnostics,  specific  ioniza¬ 
tion  stages  (or  ionic  species)  need  to  be  considered.  A 
scheme  suitable  for  application  in  a  wide  range  of  plasma 
conditions  is  examined.  Connections  are  made  to  the  non¬ 
linear  screening  model  mentioned  previously.  Implementa¬ 
tion  of  a  computer  code  incorporating  the  above  algorithm 
is  currently  in  progress.  This  should  provide  DNA  with  a 
suitable  data  base  for  a  variety  of  atomic  calculations  in 
dense,  hot  plasmas. 
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SECTION  2 

THE  EQUATION-OF-STATE  FOR  DENSE,  HOT  PLASMAS 

In  developing  a  theoretical  model  for  the  computation 
of  equation-of  state  (EOS)  data  and  other  thermodynamic 
properties  of  plasmas  over  a  wide  range  of  densities  and 
temperatures,  one  is  faced  with  the  task  of  properly 
describing  the  physics  of  vastly  different  system  condi¬ 
tions  and  properties.  In  the  low  density,  high  temperature 
region,  the  plasma  is  "weakly  coupled" — that  ,  he  random 

thermal  motion  dominates  over  the  Coulomb  :  enaction 
between  charged  plasma  particles.  In  this  rrgi  n,  many- 
body  field  theoretic  methods  can  be  used  to  obtain  the 
EOS-data.  However,  for  high  densities  and  low  temperatures 
(as  obtained  for  shock-compressed  materials,  the  superdense 
core  of  laser-fusion  plasmas,  laser-imploded  solid  targets, 
etc.),  the  plasma  is  in  the  "strongly-coupled"  condition. 
Coulomb  interaction  dominates  over  the  thermal  motion  of 
the  plasma  particles  and  perturbative  methods  are  no 
longer  applicable. 

In  this  context,  the  local  density  functional  scheme 
3 

of  Kohn  and  Sham  ,  generalized  to  arbitrary  temperatures 

2 

by  us  ,  offers  a  very  useful  and  practical  method  of 
calculation  of  EOS-data.  We  have  developed  such  a  scheme 
for  two-component  systems  with  arbitrary  coupling  condi¬ 
tions  applicable  to  both  pure  plasmas  (partially  ionized 
Aluminum,  for  example)  as  well  as  for  plasmas  seeded  with 
impurity  ions. 

In  this  density  functional  method  (DFT) ,  the  charge- 
density  profiles  of  ions  and  electrons  of  the  inhomogeneous 
plasma  play  a  central  role.  In  a  partially  ionized  plasma 
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at  an  arbitrary  density  and  temperature,  the  positive 
ion-centers  form  an  attractive  potential  for  the  electrons 
some  forming  bound  states,  others  occupying  continuum 
states.  In  the  DFT,  both  the  bound  and  continuum  states 
are  treated  quantum-mechanically  via  solution  of  the 

II 

Schrodinger  equation.  This,  in  turn,  accounts  for  a 
significant  improvement  over  other  semi-classical  approxi¬ 
mations  commonly  used  for  plasma  calculations'*  .  The 
effective  potential  in  the  DFT  contains,  aside  from  the 
nuclear  and  electrostatic  terms,  an  exchange-correlation 
contribution  arising  from  Coulomb  interactions  between 
charged  particles.  The  appropriate  degeneracy  effects 
are  maintained  through  Fermi-Dirac  statistics  for  the 
electrons.  The  electronic  charge-densities  are  con¬ 
structed  from  the  probability  densities  obtained  from  the 
wave  functions  and  are  weighted  by  the  occupation 
probabilities  of  the  various  electronic  states  according 
to  F-D  statistics.  The  effective  potential  is  then 
constructed  from  the  charge  density  profiles  (with  its 
associated  exchange-correlation  potential)  and  the 
procedure  iterated  to  self-consistency.  The  self- 
consistent  electron  and  ion  charge  density  profiles  are 
used  to  calculate  various  thermodynamic  properties  and 
EOS-data  (pressure,  total  energy,  specific  heat,  etc.). 

The  mathematical  details  of  the  scheme  are  given  in 
Appendix  I. 

At  this  point,  the  following  aspects  of  the  scheme 
should  be  noted: 

1.  The  well-known  inadequacies  of  the  Thomas- 
1  4 

Fermi  method  '  ,  extensively  applied  for  EOS-calculations , 
are  removed.  That  is,  the  kinetic  energy  is  treated 
exactly  and,  unlike  the  Thomas-Fermi  method,  the  exchange- 
correlation  effects  are  included. 


2.  The  exchange-correlation  potential  V 

A  v 

needed  in  this  scheme  can  be  calculated  by  many-body 
theory  techniques.  This  has  been  developed  and  the 
relevant  Feynman  diagrams  computed  numerically^'^. 

These  quantities  are  thus  readily  available  for  use 
in  the  above  algorithm. 

3.  Through  the  effective  potential,  the  scheme  also 
provides  the  correct  dielectric  function  for  the  plasma. 

In  strongly-coupled  plasmas,  the  usual  mean-field  approxi¬ 
mations  are  no  longer  valid.  Thus,  complicated  nonlinear 
methods  are  needed  to  calculate  the  plasma  dielectric 
function.  Moreover,  the  present  scheme  is  equally  valid 
for  strongly  as  well  as  weakly-coupled  plasmas.  The 
dielectric  functions  obtained  through  this  scheme  can  be 
used  for  a  variety  of  plasma  calculations  (energy  loss  of 
fast-charged  particles  moving  through  the  plasma,  pair 
distribution  functions  and  structure  factors,  bulk  thermo¬ 
dynamic  properties,  etc.). 

In  the  "weakly  coupled"  domain,  it  is  simpler  to  use 
many-body  field  theoretic  methods  to  calculate  EOS-data. 

The  ionic  and  the  electronic  contributions  can  be 
separated  as  a  first  approximation  and  evaluated  individ¬ 
ually.  Except  for  very  low  temperatures,  the  electronic 
contributions  dominate  the  overall  EOS-properties .  At  a 
given  density  and  temperature,  the  interacting  electron 
system  pressure  (or  total  energy)  can  be  calculated  by 
evaluating  the  kinetic  energy  part  exactly  and  the  inter¬ 
action  energy  part  through  a  many-body  perturbation  method. 
The  first  order  exchange  term,  arising  from  the  Pauli 
exclusion  principle,  is  purely  quantum-mechanical  in  nature 
The  correlation  energy  is  evaluated  in  this  scheme  via 
numerical  summation  of  the  relevant  Feynman  (ring)  diagrams 
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The  "ring"  diagram  sum  is  known  to  yield  the  Gell-Mann 

g 

and  Brueckner  result  in  the  degenerate  limit  and  goes 

9 

over  to  the  Debye  ,  '  lue  for  the  correlation  energy  in 
the  non-degenerate  (classical)  limit.  Thus  they  are 
expected  to  interpolate  properly  through  the  intermediate 
degeneracy  region.  The  mathematical  details  are  given  in 
Appendix  I . 

The  numerical  results  for  the  electronic  pressure 

24 

obtained  from  this  scheme  for  electron  densities  10 
20  -3 

10  cm  and  temperatures  0  -  10  eV  are  given  in  Table  1. 
The  drastic  difference  from  the  ideal  gas  EOS  are  illus¬ 
trated  in  Figure  1.  The  differences  progressively  reduce 

with  decreasing  density  and  increasing  temperature,  but 

20  -3 

even  for  a  density  as  low  as  10  cm  and  a  temperature 
of  1  eV  (=10,000°K),  the  difference  is  still  more  than 
50%.  This  shows  that  the  ideal  gas  EOS  is  not  valid  in 
this  density-temperature  range.  At  high  density  and  low 
temperatures  the  exchange  effect  dominates  the  interaction 
energy,  whereas  in  the  moderate  density  and  temperature 
ranges,  the  correlation  energy  is  more  important.  As  the 
temperature  increases,  the  very  slow  rise  of  the  pressure- 
curve  for  a  fixed  density  is  because  the  increase  in  the 
correlation  contribution  is  partly  offset  by  the  decrease 
in  the  exchange  effect. 

These  results  indicate  that  there  exists  a  wide 
range  of  densities  and  temperatures  where  systematic 
investigations  of  the  EOS-properties  are  needed  to  avoid 
drastic  errors  in  the  data.  Even  in  the  weakly-coupled 
domain,  many-particle  effects  are  very  substantial.  For 
strongly-coupled  plasmas,  ion-ion  and  ion-electron  inter¬ 
actions  are  important  in  addition  to  electronic  contribu¬ 
tions.  Partially  ionized  plasma  conditions  present  further 


complications.  The  model  outlined  here  incorporates  these 
effects  in  a  systematic  way.  The  computer  code  being 
developed  based  on  this  algorithm  is  expected  to  provide 
accurate  EOS-data  over  a  wide  range  of  plasma  conditions. 


SECTION  3 


THE  PLASMA  POLARIZATION  EFFECT 


The  screening  of  a  static  external  charge  immersed 
in  a  plasma  has  been  of  central  importance  not  only  in 
plasma  physics  but  also  in  astrophysics  and  semi-conductor 
physics.  Traditionally  the  static  screening  has  been 
described  in  terms  of  a  shielded  potential  of  the  form 
V  (r)  =  —  e  7  (  >  is  the  screening  parameter) .  Debye 

*3  IT 

Huckel  and  Thomas-Fermi-type  models  fall  under  this 
category.  Physically  the  electric  field  of  the  external 
charge  polarizes  the  surrounding  medium  so  that  unlike 
charges  pile  up  around  the  external  charge  whereas  like 
charges  are  repelled  from  it.  This  results  in  modifying 
the  bare  Coulomb  potential  to  a  "screened"  potential.  A 
simple  approximate  form  of  the  screened  potential  is  given 
by  Vs. 

A  general  approach  to  the  plasma  polarization  effect 
can  be  made  via  an  appropriate  dielectric  response  function 
In  general,  this  approach  leads  to  complicated  nonlinear 
methods'*-0  which  are  very  difficult  to  apply  to  plasma 
problems.  Within  linear  response  theory,  however,  the 
random  phase  (RPA)  dielectric  function  provides  a  useful, 
practical  approach. 

In  recent  years,  the  RPA  dielectric  function  was 
generalized  to  arbitrary  temperatures  and  densities  and 
a  more  complete  analysis  of  static  screening  was  done'*-'*'. 
This  analysis  showed  that  the  effective  screened  potential 
goes  over  to  the  Debye  potential  in  the  high  temperature, 
low  density  region  and  to  the  Thomas-Fermi  potential  in 
the  low  temperature,  high  density  region.  The  RPA-screened 


potential  thus  provides  a  natural  way  to  interpolate 
between  these  two  limits  and  provides  a  scheme  to  investi¬ 
gate  the  screening  effects  in  the  intermediate  range  of 
densities  and  temperatures.  Here  neither  the  Debye  nor 
T-F  model  is  applicable.  If  the  temperature  of  the  plasma 
is  of  the  order  of  the  Fermi  temperature,  T„,  drastic 

r 

differences  were  obtained  in  the  level  spectrum  calculated 
using  the  RPA-screened  potential  as  compared  to  the  spec¬ 
trum  obtained  from  the  Debye  potential.  Degeneracy  effects 
and  a  correct  momentum  dependence  account  for  these 
differences.  This  is  illustrated  in  Table  2. 

Within  linear  response  theory,  the  RPA-screened 
potential  describes  the  screening  effects  better  than  the 
Debye  or  T-F  models.  The  linear  response  theory,  however, 
underestimates  the  charge  pile-up  around  the  ions.  This 
produces  smaller  screening  shifts  of  the  ionic  energy 
levels  than  would  be  obtained  from  a  more  accurate  model. 

A  full  nonlinear  self-consistent  model  is  therefore 
required  for  appropriate  treatment  of  screening  effects. 

The  model  we  have  developed  treats  the  nonlinear  screening 
effects  exactly.  The  correct  charge  density  profile 
around  an  ion  is  constructed  from  the  probability  densities 
obtained  from  the  wave  functions  and  is  weighted  by  the 
occupation  probabilities  of  the  various  electronic  states 
according  to  Fermi  statistics.  The  effective  atomic 
potential  contains,  aside  from  the  nuclear  and  electro¬ 
static  terms,  an  exchange-correlation  contribution  arising 
from  the  Coulomb  interactions  between  charged  particles. 
These  last  two  quantities  depend  on  the  effective  atomic 
potential  via  solution  of  the  Schrodinger  equation.  The 
procedure  is  iterated  to  self-consistency.  Such  a 
numerical  calculation  is  approximately  ten  times  faster 
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than  a  self-consistent  Hartree-Fock  calculation.  The 
model  has  the  merit  of  treating  the  screening  effect 
exactly  and  the  many-body  effects  in  an  adequate  way. 

Our  model  was  applied  to  Neon  and  Argon  impurities 
in  a  dense,  hot  hydrogen  plasma.  At  the  chosen  temperature 
and  density  (T  =  200  eV,  ng  ^  10^  -  10^  cm  ,  hydrogen 
is  fully  ionized.  A  comparison  of  the  results  (Table  3) 
of  the  self-consistent  calculation  with  those  from  the 
Debye  model  show  that  the  energy  levels  obtained  through 
the  Debye  potential  are  drastically  deeper  than  the  self- 
consistent  spectrum.  This  is  because  the  linear  screening 
model  (Debye-Huckel)  underestimates  the  charge  pile-up 
around  the  ion.  The  correct  charge  density  profile  is 
obtained  from  the  self-consistent  potential  which,  in 
turn,  produces  larger  screening  shifts.  The  other  con¬ 
sequence  of  this  effect  is  that  the  Debye  potential 
incorrectly  predicts  a  larger  number  of  bound  states. 

Also,  the  line  shifts  obtained  in  the  Debye  model  are 
substantially  different  from  the  self-consistent  model 
results.  This  is  important  because  the  bound-bound 
transition  probability  varies  as  the  cube  of  the  level 
separation.  The  usefulness  of  our  model,  which  yields 
the  correct  screening  effect  in  dense,  hot  plasmas,  is 
therefore  obvious. 

The  existence  of  different  ionic  species  (hydro- 
genic,  helium-like,  etc.)  in  experimentally  observed 
spectra  from  hot  plasmas  is  an  additional  aspect  that 
needs  to  be  considered.  For  determination  of  the  correct 
energy-level  spectrum  of  one-,  two-  or  multi-electron 
ions  in  hot  plasmas,  the  above  model  needs  to  be  extended. 
For  purposes  of  illustration,  the  case  of  a  1-electron 
hydrogen-like  ion  immersed  in  hot  plasma  is  considered 


« 


X 

i 

i 


in  Appendix  IIB.  For  most  practical  purposes,  the  1- 
electron  spectrum  is  within  a  few  percent  of  the  level 
spectrum  obtained  in  the  self-consistent  calculation 
discussed  above. 

12 

Previous  calculations  by  Skupsky  and  by  Davis  and 

Blaha"^  utilize  two  different  approaches  to  the  problem 

of  a  1-electron  spectrum  in  dense,  hot  plasmas.  Davis  and 

Blaha"^  employ  an  approach  similar  to  a  self-consistent 

Hartree-Fock  method  and  leave  out  the  self-interaction  term 

from  the  screened  potential.  Also,  the  exchange-correlation 

potential  was  included  only  for  the  continuum  electrons  and 

not  for  the  bound  electrons.  Calculations  were  done  for 

level  shifts,  coefficients  of  transition  probabilities, 

9  + 

and  electron  collision  cross-sections  of:  Ne  for  tempera¬ 
tures  200  -  500  eV  and  electron  densities  (1-6)  x  10^  cm 

and  Ar7+  for  temperatures  1000  -  2000  eV  and  electron  den- 

25  -3 

sities  (2-8)  10  cm  .  The  general  trends  of  the  results 

indicate  substantial  screening  effects  on  those  atomic 
parameters . 

12 

In  his  scheme,  Skupsky  solves  Poisson's  equation 
self-consistently  to  obtain  an  effective  potential  which 
includes  both  bound  and  free  electrons.  However,  the  free- 
electron  density  is  evaluated  in  a  semi-classical  approxi¬ 
mation  and  many-body  effects  were  neglected.  In  contrast, 
our  model  constructs  the  free  electron  charge  density  via 
solution  of  the  Schrodinger  equation  and  includes  exchange- 
correlation  effects. 

In  conclusion,  we  note  that  our  model  treats  the 
nonlinear  screening  effect  in  an  exact  fashion  and  the 
many-body  effects  adequately.  For  plasmas  in  the  inter¬ 
mediate  degeneracy  region,  as  occurs  in  the  dense  core 
of  laser  fusion  plasmas,  shock-compressed  material  at 
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elevated  temperatures,  etc.,  this  model  can  readily  be 
applied  to  investigate  screen  effects  on  atomic  parameters 
For  multi-electron  ions  in  plasmas,  however,  this  model 
needs  to  be  extended.  Currently  we  are  in  the  process  of 
developing  such  a  model.  This  will  provide  an  important 
diagnostic  tool. 
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APPENDIX  I 


THE  EOS-MODEL 

In  this  section,  we  will  outline  our  EOS-model  and 
the  method  of  calculation.  Consider  a  partially  ionized 
plasma  (say  Aluminum  or  Argon)  at  a  given  density  and 
temperature.  Because  of  the  polarization  effect,  electrons 
will  tend  to  pile  up  around  the  ions--some  forming  bound 
states,  others  occupying  continuum  states.  Considering  a 
central  ion,  the  other  positive  ions  will  be  repelled 
from  it  whereas  electrons  will  be  attracted.  A  sphere  of 
radius  R  around  the  central  ion  will  be  electrically 
neutral,  i.e. 


Z  -  jp^(r)dr  +  Z  |Op(r)dr  =  0 


(1-1) 


where  o„(r)  and  o  (r)  are  the  electron  and  ion  density 
e  p 

profiles  around  the  central  ion.  Note  that  Eq.  ( I  —  1 )  really 

defines  the  radius  R.  If  n,  is  the  mean  number  of 

b 

electrons  bound  to  the  central  ion. 


(1-2) 


represents  the  average  charge  per  ion.  n^  is  defined  by 

integrating  the  R.H.S.  of  Eq.  (I-8g)  over  all  r.  According 

to  density  functional  theory,  the  thermodynamic  potential 

71  of  the  system  is  uniquely  determined  if  the  correct 

density  distributions  n  (r)  and  p  (r)  are  known: 

e  P 


where  T 


M 


is  the  kinetic  energy  sum  of  ions  and 


electrons,  and 


e  '  '/I 

i  ■  iff 


P  (r ) dr  + 
e 


pp(r)dr  + 


1  f  pe(?)pe(?' 

2J  |  — ’ 

Z2fPe(^)pe(?' 

2J  |F^r- 


)  drdr 1 


+  /  F  p  dr 
/  x  c  e 


-  u  IP. (r) dr  (1-4) 
e  /  e 


)  drdr ' 


-  f pe(^  n(r' 5  d^r'  f  ei 
=  -Z  /  — - 2 -  +  Fe  p  , 

;1  J  |F^?'|  J  c  e 


p  dr 
P 


-  u.  p  (r) dr  (1-5) 
1  P 


p  dr 
P 


(1-6) 


In  the  above,  F  is  the  exchange-correlation  free  energy 

X  c 

of  the  electrons  and  is  the  electron  chemical  potential. 

The  exchange  free  energy  of  ions  is  negligible.  Therefore 

only  the  correlation  part  F  is  retained.  To  simplify 

c  ei 

calculations,  we  will  neglect  Fc  which  is  negligibly 
small  in  most  cases  of  interest. 

To  take  into  account  the  charge-neutrality  condition 
(1-1),  we  consider  the  quantity 


-■[ 


Z  +  Z  J  pe(r)dr 


-I 


Pp (r)dr 


(1-7) 


where  A  is  Lagrange  multiplier.  The  stationary  condition 
of  PT  leads  to  the  following  system  of  equations  for  the 
electrons  and  ions  to  be  solved  self-consistently : 
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with 


o  ( r 
e 


,T)  =  Y.  Ui 


(I-8a) 


,BtEi  ’  ue>  +  1 


(the  sum  over  i  involves  the  sum  over  bound  and  continuum 
6  6 

states)  and  V  =  5F  /op  (r)  is  the  excxiange-correlation 

XC  X  C  6 

potential  for  the  electrons. 

For  the  ions,  we  get  a  similar  set  of  equations: 


(r' )  -  Z  c  (r ' )  dr' 


V1  (r) 
c 


Up-Z.\  «.(r)  =  £  i  ^  i  ( r ) 


(1-9) 


Due  to  their  heavier  mass,  the  ions  form  a  classical  sub¬ 
system  for  most  cases  of  interest.  Hence  the  ion  charge 
density  distribution  can  be  well  represented  by 


Op(r)  =  e 


1  exp  ^-6Vi(r) 


where  -  V  (R)  +  ZA 


and  V^(r)  =  Z 


[,  _J  k(?l>  -  Zpp(r')]dr' 


( I-9b) 


+  V1 (r)  -  V1 (R) 
c  c 


The  calculation  thus  reduces  to  obtaining  the  correct 

density  distributions  p  (r,T)  and  p  (r-T)  self- 

e  p 

consistently. 


( I  —  9  c ) 
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(I-Sb) 


Defining  u  =  u  -  V  (R)  -  \ 
J  e  e  xc 


Ea.  (1-8)  is  cast  into  the  form 


‘  X  +  Vef f (r»  °e'T) 


+J i  ( r )  =  (r) 


( I-8d) 


e  Z  I  Jpe(r’  }  ‘  zpD(r'  >]  dr' 

with  V*ff(r,  :a,T)  =  -  \  +1  -L® - £-  J 


r  -  r 1 


+  Ve  (r) 
xc 


-  (R) 

xc 


( I-3e) 


The  electron  density  distribution  is  the  sum  of 
bound  and  continuum  state  contributions: 


o  (r) 
e 


°e(r> 


oC(r)  +  o  ( R) 
e 


( I-8f ) 


where 


b  .  . 
o  (r) 
e 


Y,  2 ( 2 i + 1 ) |*i4(r) |2  f(ei,Ue) 
i  f  ?- 


( I-8g) 


CO 


( I-8h) 


o 

c ( R)  is  effectively  the  bulk  electron  density  pQ  and  the 

j's  are  spherical  Bessel  functions.  In  actual  calculations, 

it  is  convenient  to  work  with  ,  p(R)  and  the  mean 

density  o.  That  is,  we  can  start  off  with  a  given  mean 

electron  density  n  =  Z  p  and  a  chemical  potential  u 

0  0 

and  adjust  (?0  to  ensure  electrical  neutrality. 

Once  the  electron  and  ion  charge  density  distributions 

p  (r,T)  and  (r,T)  are  calculated  self-consistently , 
e  o 

all  the  thermodynamic  quantities  can  be  readily  obtained. 

For  example,  tne  pressure  P  of  the  plasma  is  given  by 
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(V  =  Volume) ,  with  i  obtained  from  Eq.  (1-3)  at  any 
density  and  temperature.  The  total  energy  E  of  the 
plasma  is  given  by 


E  =  n  +  u  N  +  U  N  +  TS 
e  e  P  P 


(1-11) 


where  the  entropy  3  is  obtained  from  the  set  of  eigen¬ 
values  .  '  s  ana  the  occupation  probabilities  f.'s  already 
computed ; 


S  =  )  f . inf.  +  (1-f . ) • n ( 1-f . ) 

<  I  1  l  l  l 

For  hydrodynamic  calculations,  several  other  quantities 
like  (dE/dp) ,  (dP/do) ,  (dP/dT)  etc.  are  needed.  These 
again,  are  readily  obtained  from  knowledge  of  the  total 
energy  E  and  pressure  P. 


(1-12) 


MANY- BODY  PERTURBATION  METHOD 


As  mentioned  in  Section  2,  the  above  scheme  is  par¬ 
ticularly  applicable  for  strongly-coupled  plasmas.  However, 
in  the  weakly-coupled  domain,  it  is  simpler  to  use  many- 
body  perturbation  theory  to  calculate  EOS-data.  This  is 
discussed  below.  Of  particular  importance  is  the  electronic 
contribtuion  to  the  EOS.  Let  us  focus  our  attention  to  an 
interacting  electron  system  of  density  at  a  temperature 

T.  The  thermodynamic  potential  P[c,T]  can  be  written  as 


[o,T]  = 


xc 


(1-13) 
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where 


is  tiio  non-interacting  part 
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dx  x  — 


,  -1  2  . 
(t  x  -a) 


(1-14) 


—2  ?  ~/3 

with  _F  =  h  /2m,  the  Fermi  energy  of  the  elec¬ 
trons,  x  =  k  /k_  .,  and  a  =  y/K^T  is  the  degeneracy 

parameter.  For  arbitrary  densities  and  temperatures,  a  is 
to  be  determined  from  the  phase  space  integral. 


(t 


r  ,  -1  2  ,  ] 

e(t  x  "a)+l 


T/T„  .  )  .  At  zero  temperature 
Fermi  r 


n 

“o  ( 
V 


(1-15) 


( I-14a) 


with  the  corresponding 

P  (t=0)  =  | os 
o  5 

At  high  temperatures, 
equation  of  state 


pressure 

F 

one  recovers  the  ideal 


gas 


P  - (I-14b) 

°  m  B 

rT1  —  .*3 

The  interaction  part  of  the  thermodynamic  potential 
dxc[o,T]  is  a  difficult  many-body  quantity.  The  exchange 
part  arises  from  tne  Pauli  exclusion  principle  whereas  the 
correlation  effects  arise  from  th,e  Coulomb  repulsion  between 
electrons  which  tend  to  keep  them  apart.  The  leading 
exchange  contribution  ft  is 
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-Vo 


I  t  X  -  c  )  .  t-1 

e  +1  e  y 


-1  2 -:t) 


At  zero  temperature,  we  have 


(1-16) 


n  (t=o> 

x 


,,  2  ,4/3 
( 3  tt  a) 


(I-16a) 


whereas,  in  the  high  temperature  limit,  T-*-<» 


V"~°>  e2  2. ,4/3  4eF 

— V - ^  <3’  »)  9KJ 


(I-16b) 


Since  in  the  weakly-coupled  regime,  the  interaction  is 
small  compared  to  the  average  kinetic  energy,  diagram¬ 
matic  perturbation  theory  can  be  used  to  calculate  the 
correlation  energy.  The  "ring"  diagram  sum  gives  the 
correlation  contribution: 


ri  's* 


+  ^X(q,vn) 

q  n 


(1-17) 


where  the  polarization  function  is 


x(q/vn)  =  -2 


f  (p+q)  -  f(p) 


( 2tt )  3  ihv  -  (e  -  e  ) 
n  p+q  p 


(I-17a) 


with  vn  =  2nTTl<BT/h,  (n  =  integers).  f  ( o )  is  the  Fermi 
function.  In  cha  high  temperature  limit,  most  of  the 


•r'f' 


0)  of 


i.? 


contribution  comes  from  the  static  part  (v  = 

Eq.  (1-17)  and  from  the  small  momentum  transfers  (q-*0)  . 
In  this  T-*-00  limit. 
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3  3/2 
e  p 


(kbt) 


1/2 


( I-17b) 


which  is  the  usual  Debye  contribution  to  the  correlation 
energy.  At  the  opposite  limit  of  zero  temperature  and 
very  high  densities,  one  obtains  the  well-known  Gell-Mann 

g 

and  Brueckner  result. 
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APPENDIX  IIA 


THE  PLASMA  POLARIZATION  EFFECT 


Let  us  consider  a  homogeneous  plasma  of  density  n 
at  a  temperature  T.  If  an  external  impurity  charge  +Ze 
is  immersed  in  the  plasma,  it  will  polarize  the  medium 
and  bring  about  a  redistribution  of  charge  density.  Within 
linear  response  theory,  the  effective  potential  V  (r) 
at  a  distance  r  from  the  static  impurity  charge  +Ze 
is  given  by 

2  00 

Veff(r,n,T)  =  J  dq 

o 

q 


Sin (qr ) 


qr 


1  - 


47ie‘ 


(II-l) 


(q ,  n ,  T) 


where  ;<(q/n/T)  is  an  appropriate  density  and  temperature 
dependent  response  function.  The  random  phase  form  (RPA) 
of  x (q*n,T)  generalized  to  arbitrary  temperature  gives 


(II-2) 


where  q  and  k  are  momentum  variables,  f (k)  is  the  Fermi 

2  1/3 

function,  and  k^  is  the  Fermi  momentum;  k_  =  (3tt  n)  . 

±  r 

Veff  satisfies  the  consistency  condition  that  the  total 
induced  charge  density  must  compensate  the  external  charge, 


( 1 1-3 ) 


Unlike  the  Thorna.s-Fermi  theory,  this  induced  charge  density 
is  finite  everywhere,  including  the  origin. 

One  can  show  that  the  q-*-0  limit  of  Eq.  ( 1 1  —  1 )  is 
the  static  screened  Coulomb  potential  (SSCP) : 


Veff(r'n-T)  5^ 


Ze  -r5(n,T) 
r 


( I 1-4 ) 


where  the  parameter  t,  is  defined  by 

2  ?  7  2  4kF 

S2(n,T)  =  C"’TF  /  dk  f(k);  -TF  =  ~ 

*  o 

o 

Further,  as  T-*°°,  Sq.  (II-l)  yields  the  Debye-Huckel 
potential 


Veff(r’a'T)  ^5- 


Ze2  -r£;_  ,2  4TTne2 

—  e  D'  -d  =  ITT" 

B 


(II-5) 


( II-6 ) 


In  the  opposite  limit  of  T  =  0,  one  obtains  the  T-F 
screened  potential 


Vef£(r'n'T) 


— —  e“r£ TF 
r 


( 1 1  —  7 ) 


Thus,  the  Deb'-'e  and  the  Thomas-Fermi  screening  lengths 
are  seen  to  be  limiting  forms  of  the  generalized  screening 
parameter  £(n,T)  which  contains  the  appropriate  degeneracy 
effects.  The  full  momentum  (q)  dependence  of  veff^r^  in 
Eq.  (II-l)  yields  a  further  difference  from  the  Debye  or  the 
T-F  potential. 

Level  shifts  of  bound  electronic  energy  levels, 
supported  by  an  attractive  impurity  potential  in  a  plasma, 
has  been  investigated  extensively.  The  radiation  from 
electronic  transitions  between  these  shifted  energy  levels 
is  used  as  a  plasma  diagnostic.  The  problem  requires  the 
solution  of  the  Schrodinger  equation 


1.1 

l  2m 


+  +  v  (r 
,  2  _  2  eff' 

dr  2mr 


,n,T!)  Rn 


£(r>  ’  V  V(r)  (II-8> 
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where 


(r)  represents  the  radial  part  of  the  electron 


wave  function  with  energy  eigenvalue  e  Many  previous 

calculations  have  used  the  Debye  potential  in  Eq.  (II-8) . 
Our  results  utilizing  the  R.P.A.  screened  potential  in 
Eq.  ( 1 1— 8 )  and  obtaining  the  energy  levels  by  numerical 
solution  of  the  Schrodinger  equation  is  given  in  Table  2. 

The  linear  screening  method,  however,  underestimates 
the  charge  pile-up  around  a  given  ion.  A  full  treatment 
therefore  requires  a  nonlinear  self-consistent  method  to 
investigate  plasma  screening  effects.  We  will  illustrate 
our  method  for  tne  general  case  of  a  partially  ionized 
plasma  at  an  arbitrary  temperature  T. 

If  p  (r,T)  and  p  (r,T)  represent  the  appropriate 
e  p 

electron  and  ion  density  profiles,  then  the  requirement  of 
self-consistency  implies  the  following  set  of  coupled 
equations : 


Xr-  -  -  + 

2  r 


p  (r'  , T)  -  Zp  (r '  ,T)1  dr' 

— - E - - -  +  V  (r,T)  i p.  (r) 

I r  -  r'  I  xc 


=  e  (r ) 


with  oe(r,T)  =  Y.  I»i(r) I2  — 

i  e^  1  +1 


( II-9 ) 


(11-10) 


(the  sum  over  i  implies  the  sum  over  bound  and  continuum 
states) ,  and , 

72  77  fkXr'.V  ~  Zp  (r  ’  ,T)1  dr '  . 

-_  +  _.  ij  J^= - - J -  +  vc(r)Jn(r) 


=  e  ■  cf  .  (r ) 
li 


(11-11) 
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In  the  above,  Ve  (r,T)  is  the  exchange-correlation 

xc  ^ 

potential  for  electrons,  V  (r)  is  the  correlation 
potential,  and  Z  is  the  effective  charge  of  a  partially 
ionized  ion.  In  most  practical  cases  of  interest  the  ions 
can  be  treated  classically  because  of  their  heavier  mass. 

Thus  the  solution  of  Eq.  (II-l)  becomes: 

-6V. (r) 

o  (r,T)  =  o  e  1  (II-12) 

P  ° 


In  setting  up  the  above  equation  many-body  effects 
have  been  included  via  Ve  (r,T)  for  the  electrons  and 
V1  for  the  ions.  Note  that  the  induced  charge  displace¬ 
ment  does  compensate  the  nuclear  charge  over  a  large 
distance,  thus  accounting  for  complete  screening.  The 
above  method  was  applied  to  Neon  and  Argon  impurities  in 
a  fully  ionized  hydrogen  plasma.  The  results  are  shown 
in  Tables  3(a)  and  3(b). 
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APPENDIX  IIB 


IIYDROGENIC  IONS  IN  PLASMAS 


For  purpose  of  illustration,  we  consider  the  specific 
case  of  a  hydrogen--like  ion  (i.e.,  an  ion  with  a  single 
electron  bound  to  it)  in  a  plasma.  (The  extension  to 
multi-electron  ions  is  straightforward  but  the  details  will 
not  be  given  here) .  A  rigorous  approach  would  require  the 
self-consistent  solution  of  a  non-local  integro-dif f erential 
equation  of  the  following  form: 


(r  ,r  1  ,E)(p  (r*) 


4>i  (r) 


=  e  .  .  (r) 

11 


where  2_,  (r,r',E)  is  the  non-local,  energy-dependent  "mass" 

operator  which  involves  knowledge  of  {  (r)>  ,  the 

energy-levels,  and  the  orbitals  of  interest.  The  following 
scheme  substantially  reduces  the  complexity  of  such  a 
calculation  while  still  retaining  most  of  the  many-body 
effects . 

Let  0  n (r)  be  the  bound-state  wave  function  of  the 
n 

electrons  bound  to  the  ion  Z.  (r)  could  be  any  of  the 

Is,  2s,  2p  ...  sequence  of  states.  Based  on  physical 
considerations,  we  invoke  the  argument  that  the  screening 
term  affecting  the  bound  electron  does  not  include  the 
part  associated  with  the  bound  electron  itself.  In  other 
words,  we  exclude  the  self-interaction  effect.  The  energy 
level  spectrum  is  then  determined  by  the  solution  of 


z  rLp(r  >T)  -  i*ng(r 

2  r  J  \r  -  ?'| 


)  dr' 


jp(r,T) 


-  Vxc  l*nl(r)  1  'T  V(r)  =  £n^n£(r) 


Note  that  one  has  a  different  equation  for  each  of  the 
n£  states  (Is,  2s,  sp  ...  etc.).  The  advantage  of  this 
scheme  is  that  p(r,T)  is  the  self-consistent  charge 
density  profile  already  obtained  from  previous  nonlinear 
screening  calculations  for  the  same  density  and  tempera¬ 
ture.  Note  that  the  continuum  charge  density  pc(r,T) 


satisfies  the  requirement 


Pc(r,T)dr  =  (Z  -  1) 


Since  p(r,T)  =  Pc(r,T)  +  |<j>n£(r)| 

and  all  $  (r)  are  normalized  to  unity: 


<J>  (r)  I  dr  =  1 


TABLE  1 


Electronic  Contribution  to  Equation  of 
State  Pressure  in  Ryd/cm3 

(The  difference  from  ideal  gas  EOS  shown  in  parenthesis) 

Temperature  (eV) - *• 

Electron 

density 


(cm  3  ) 

0 

0.0316 

0.316 

3 .16 

10 

io24 

t  2.003E24 

1  (2. 003E24) 

2 . 004E24 
(2. 002E24) 

2.006E24 
(1 . 983E24) 

2.025E24 

(1.793E24) 

2 . 275E24 
(1 . 540E24) 

10  23 

j  0. 7074E23 
)  (0. 7074E23) 

0. 7076E23 
(0. 7053E23) 

0. 7085E23 
(0.6853E23) 

0. 8225E23 
(0. 590E23) 

1 . 142E23 
(0 .407E23) 

1022 

1  0.297E22 
1(0. 297E22 

0. 299E22 
(0. 296E22 

0. 302E22 
(0.279E22) 

0. 443E22 
(0.211E22) 

0. 889E22 
(0.154E22) 

1'J21 

|  C . 140E21 
IOJ.  14  0E21) 

0. 155E21 
(0.153E21) 

0. 175E21 
(0. 151E21) 

0.321E21 
(0. 089E21) 

0. 784E21 
(0. 049E21) 

1  0  .  C  6  7  E  2  0 
l  0  .  0 1  7  K  2  0 ) 

0. 078E20 
(0. 076E20) 

0.095E20 
(0. 071E20) 

0. 259E20 
(0. 027E20) 

0 . 750E20 
(0. 016E20) 

TABLE  2 


Is,  2s,  2p  Energy  Levels  for  a  Neon  Impurity  in  a 
Hydrogen  Plasma  at  T  =  100  eV  for  Three  Different 
Electron  Densities. 

(The  numbers  in  parenthesis  show  the  percentage 
differences  from  corresponding  Debye  results. 
Units  in  Rydbergs.) 


Electron 

Potential 

Bound  States 

Density 

Used  in 

Line 

(cm‘^) 

Eq. ( II-3) 

Is 

2s 

2p 

Shift 

VD 

Vf 

-72.816 

-5.316 

-4.076 

6.260 

4. 54xl024 

-75.499 

-6.623 

-5.560 

5.061 

(3.7%) 

(24.6%) 

(36.4%) 

Vef  f 

-78.695 

-7.383 

-6.591 

2.896 

(8.1%) 

(38.9%) 

(61.7%) 

24 

2.47x10 

VD 

vr 

-79.369 

-80.605 

-8.747 

-9.487 

-7.939 

-8.759 

3.570 

3.154 

s. 

(1.6%) 

(8.5%) 

(10.3%) 

Veff 

-82.800 

-9.930 

-9.330 

1.530 

(4.3%) 

(13.5%) 

(17.5%) 

Q  y* 
>  > 

-86.482 

-13.406 

-13.022 

1.540 

24 

1.0x10 

-86.842 

-13.668 

-13.302 

1.460 

s 

(0.4%) 

(2%) 

(2.2%) 

V  ** 
ef  f 

-88.043 

-14.037 

-13.781 

0.738 

(1.8%) 

(4.7%) 

(5.8%) 
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Energy-Level  Spectrum  (Atomic  Unit)  for  a  Neon  Impurity 
in  a  Hydrogen  Plasma  at  T  =  200  eV  for  Three  Different 
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Spectrum  for  an  Argon  Impurity  Ion  in  a  Hydrogen 
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CD 
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illustrated . 


t=  T/ 


Fermi 


Total  electronic  exchange-correlation  contribution 
to  the  thermodynamic  potential  ft  as  a  function  of 
temperature  for  two  different  electronic  densities 
The  dotted  curve  indicates  the  rapid  decay  of  the 
exchange  effect  with  increasing  temperature.  The 
high  temperature  part  of  ft^c  is  completely  domi¬ 
nated  by  correlation  effects. 


200  eV 


Figure  4.  The  strong  density  dependence  of  the  self-consistent  a t om i 
potential  F,q.  ( T I  —  9 )  for  a  Neon  ion  embedded  in  a  hydrogen 
plasma  at  T  -  200  oV. 


2  00  eV 

i  n24  -3 

10  cm 
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Figure  5.  Energy-level  spectrums  for  a  Neon  ion  m  a 
hydrogen  plasma.  Self-consistent  results 
using  VQff  yield  drastic  differences  from  the 

corresponding  Debye  spectrum.  This  emphasizes 
the  importance  of  the  nonlinear  screening 
method  used  in  our  model. 
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